International Journal of Applied Power Engineering (IJAPE) 

Vol. 9, No. 1, April 2020, pp. 58-66 

ISSN: 2252-8792, DOI: 10.11591/ijape.v9.il.pp58-66 


58 


□ 

Transient development of MHD natural convection flow in 
vertical concentric annulus 


Basant K. Jha, Taiwo S. Yusuf 

Department of Mathematics, Ahmadu Bello University, Nigeria 


Article Info 


ABSTRACT 


Article history: 

Received Mar 29, 2019 
Revised Apr 20, 2019 
Accepted Mar 12, 2020 


Keywords: 

Constant heating 
Hartmann number 
Implicit finite difference 
Natural convection 
Riemann-sum approximation 


This paper examines the role of magnetic field on fully developed natural 
convection flow in an annulus due to symmetric of surfaces. The transport 
equations concerned with the model under consideration are rendered 
non-dimensional and transformed into the ordinary differential equation 
using Laplace transform technique. The solution obtained is then transformed 
to time domain using the Riemann-sum approximation approach. The governing 
equations are also solved using implicit finite difference method so as to 
establish the accuracy of the Riemann-sum approximation approach at 
transient as well as at steady state solution. The solutions obtained are 
graphically represented and the effects of pertinent parameters on the flow 
formation are investigated in detail. The Hartmann number ( M ), is seen to 
have a retarding effect on the velocity, skin-frictions and the mass flow rate. 
Also, skin-friction at both surfaces and the mass flow rate within the annulus 
are found to be directly proportional to the radii ratio (2). 
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NOMENCLATURE 

B 0 Constant magnetic flux density 

C p Specific heat at constant pressure 

g Gravitational acceleration 

k Thermal conductivity of the fluid 

M Hartmann number 

Pr Prandtl number 

Q Dimensionless mass flow rate 

7T Radius of the inner cylinder 

r 2 Radius of the outer cylinder 

r' Dimensional radial coordinate 

R Dimensionless radial coordinate 

t Dimensionless time 

t' Dimensional time 

T 0 Reference temperature 

T w Temperature at both surfaces 

u' Axial velocity 

U Dimensionless axial velocity 
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Greek letters 

9 Dimensionless temperature 

v Fluid kinematic viscosity 

r Skin-friction 

a Electrical conductivity of the fluid 

p Density 

A Radii ratio 

/? Coefficient of thermal expansion 

1. INTRODUCTION 

In recent times, the natural convection with simultaneous heat and mass transfer flow of 
magnetohydrodynamics fluid has attracted many researchers in view of its various applications in 
astrophysics, geophysics, meteorology, aerodynamics, magnetohydrodynamics power generators and pumps, 
boundary layer control energy generators, accelerators, aerodynamics heating, polymer technology, 
petroleum industry, purification of crude oil, and in material processing such as metal forming, continuous 
casting wire, magnetohydrodynamics (MHD) bearings, MHD flow meters, geothermal energy excitations 
and plasma controls. Others include the study of stellar and solar structures, interstellar matter, and radio 
propagation through the ionosphere and so forth [1-9], 

The study of MHD flow problem in an annulus was first discussed by Globe [10] who considered 
fully developed laminar MHD flow problem in an annulus. In his work, he studied the problem of steady 
flow of an electrical conducting incompressible fluid in an annular space between two infinitely long circular 
cylinders under the radial impressed magnetic fluid. Later on, Georgantopoulos and Goudas [11] studied 
the free convection effects on the hydromagnetic oscillatory flow in the Stokes problem past an infinite 
porous vertical limiting surface with constant suction. Again, the analytical solutions for transient fully 
developed natural convection in open-ended vertical concentric annulus was presented by Al-Nimr [12], 
Later on, Al-Nimr [13] carried out analytical solutions for fully developed MHD natural-convection flow in 
the open-ended vertical concentric porous annulus. Al-Nimr and Database [14] presented the closed forms on 
transient fully developed free convection solutions, corresponding to four fundamental thermal boundary 
conditions in the vertical concentric annulus, while Sheikholeslami and Gorji-Bandpy [15] examined 
the numerical solution for free convection of ferrofluid in a cavity heated from below in the presence of 
external magnetic field. Sheikholeslami et al. [16] investigated MHD natural convection of nanofluid in 
a concentric annulus between a cold outer square cylinder and a heated inner circular. 

In other work, Singh et al. [17] studied natural convection in vertical concentric annulus under 
a radial magnetic field, where they observed that both velocity and temperature are more in case of 
isothermal condition compared with constant heat flux case when gap between cylinders is less or equal to 
radius of inner cylinder. Nirmal et al. [18] studied an exact solution for unsteady magnetohydrodynamic free 
convection flow with constant heat flux and concluded that the magnetic field has a retarding effect on 
the velocity while the skin-friction at the plate increases with it. Furthermore, Jha et al. [19] considered fully 
developed MHD natural convection flow in a vertical microchannel with the effect of the transverse magnetic 
field in the presence of velocity slip and temperature jump at the annular micro-channel, they established that 
increase in curvature radius leads to an increase in the mass flow rate. 

In a related article, Abbas et al. [20] investigated application of drug delivery in magnetohydrodynamics 
peristaltic blood flow of nanofluid in a non-uniform channel, while Bhatti et al. [21] considered combined 
effect of magnetohydrodynamics and partial slip on peristaltic blood flow of Ree-Eyring with wall properties. 
An analysis to investigate the combined effects of heat and mass transfer on free convection unsteady 
magnetohydrodynamics (MHD) flow of the viscous fluid embedded in a porous medium was later presented 
by Ali et al. [22], Most recently, two-dimensional magnetohydrodynamic flow of a viscous fluid over 
a constant wedge immersed in a porous medium was studied by Kudenatti et al. [23], they found that MHD 
effects on the boundary layer are exactly the same as the porous medium in which both reduce the boundary 
layer thickness. This paper is devoted to investigate the role of magnetic field on fully developed natural 
convection flow of an incompressible and electrically conducting fluid filled between two vertical coaxial 
cylinders, when the two cylinders are subjected to constant heating. 


2. RESEARCH METHOD 

Consider a transient laminar fully developed natural convection flow of an incompressible, viscous 
and electrically conducting fluid in an annulus of infinite length under the influence of transverse magnetic 
field. The z-axis is taken along the axis of the cylinder in the vertically upward direction and r'-axis is in 
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the radial direction. A magnetic field of strength B 0 is assumed to be uniformly applied in the direction 
perpendicular to the direction of flow. In the present physical situation, the inlet fluid temperature is 
maintained at T 0 , while a constant uniform heating of T w is applied at the outer surface of the inner cylinder 
and at the inner surface of the outer cylinder such that T w > T 0 as presented in Figure 1. The flow is assumed 
to be fully developed both thermally and hydrodynamic ally, and the viscous dissipation, radiation, 
and compressibility effects are neglected. Following the work of Iha et al. [24] the momentum and energy 
equations governing the present physical situation is given by (1) and (2). 



Figure 1. Schematic diagram of the problem 


du' d 2 u' 1 du' 
dt' V dr’ 2 + r' dr' 


oB 2 


P 


u'+gPO"-T, o) 


dT' 

~dt 7 


k 

PCp 


d 2 T 1 dT' 
dr' 2 + r' dr' 


The relevant dimensional initial and boundary conditions are; 


t' < 0 u' — 0, T' = T 0 for r x < r' < r 2 
1-' f u ' = °- r = T w at r' = r x 

c - u 1 u' = 0, T' = T w at r' = r 2 


Introducing the following dimensionless quantities in (1) and (2). 


t = ^ = = ^ , M 2 

U =u'v[g/3(.T w 1 -Wrl 2 ]- 1 


aB 0 2 " r l 2 g _ ( T '- T 0 ) p r _ rfv_ 

pv ’ (r w -r 0 ) ’ k 


Equations (1) and (2) in dimensionless form are obtained as follows: 


du 

at 


= Pl + -T-~ M2 U + 9 

dR 2 R dR 


„ ae a 2 e , i ae 
Pr — =-1- 

at dR 2 R dR 


The initial and boundary conditions in dimensionless form are: 
t < 0 U = 0, 0 = 0 for 1 < R < A 


( 1 ) 

( 2 ) 


(3) 


(4) 


(5) 

( 6 ) 

(7) 
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t > 0 


r U = 0, 
1 U = 0, 


U = 0, 0 = 1 
0 = 1 


at ft = 1 
at ft = 2. 


( 8 ) 


The physical quantities used in (1) to (6) are defined in the nomenclature. The solution of (5) and (6) 
with the associated initial and boundary conditions (7) and (8) can be obtained by using the Laplace 
transform technique. Defining the following transform variables. 


t/(ft,s) = f™U(R,t)e st dt, 6(R,s) — / o °°0(ft,t)e st dt, 


(9) 


Where the Laplace parameter, (s > 0) in (5) and (6) are transformed into the Laplace domain using 
the initial condition (7) to obtain. 


d 2 U , 1 dU ,,.n , n ^ 

— +-(M 2 +s)U = -0 

dR 2 R dR v ' 


d 2 8 1 dd s 

—-1- sPrd = 0 

dR 2 RdR 


Applying Laplace transform technique (9) on the boundary conditions (8), we have: 


17 = 0 , 0 = Vs 
£7 = 0, 0 = Vs 


at ft = 1 
at ft = A 


The solution of (10) and (11) in Laplace domain subject to the boundary conditions (12) are: 
U(R,s) = C 3 I 0 (RS) + CMRS) - 
0(ft,s) = C^RyfsPr) + C 2 K 0 (Ry[sPr) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


Equation (13) and (14) are to be inverted in order to obtain their solutions in the time domain. 
Due to the complex nature of these inversions, we adopt a numerical procedure used in Iha and Yusuf [25] as 
well as Jha and Apere [26] which is based on the Riemann-sum approximation. According to this technique, 
any function in the Laplace domain can be inverted to the time domain as follows: 


ft(ft, t) = f gft(ft, £ ) + Re S7=i P (ft-£ + iJ f) C-1) J ] , 1 < « < A 


(15) 


where Re refers to the real part of i = V—1 the imaginary number. N is the number of terms used in 
the Riemann-sum approximation and £ is the real part of the Bromwich contour that is used in inverting 
Laplace transforms. The Riemann-sum approximation for the Laplace inversion involves a single summation 
for the numerical process its accuracy depends on the value of £ and the truncation error dictated by M. 
According to Tzou [27], the value of et that best satisfied the result is 4.7. 

2.1. Skin friction and mass flow rate 

The skin friction at ft = 1, f^ft, s), and ft = 71, f^(ft,s)in Laplace domain is obtained by 
differentiating equations (13) and (14), respectively. While, the mass flow rate of the fluid through 

the annular gap in Laplace domain Q (ft, s), is obtained by evaluating the integral 271 ft U (ft, s)dR. 
The solutions are as follows: 


Tt 


dU I 


= 5(C 4 K^AS) - C 3 h 


s(Pr-l)-M 2 I 


- rc 1 / 1 (VsPr)-C 2 K’ 1 (VsPr)j 

L s(Pr-l)-M 2 J 

(16) 

rc 1 t 1 (AV3T)-c 2 Ki(AVST)j 

L s(Pr-l)-M 2 J 

(17) 

hdS)) -^(AK.dAS) - ^(8))]- 
^(ak^aVIpP) - K^VIP?))}} 

(18) 
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where I 0 , K 0 , I t , K 1 are the modified Bessel function of first and second kind of order 0 and 1 respectively. 
In the same manner, the solutions are inverted to the time domain by applying the Riemann-sum 
approximation stated in (15). 

2.2. Validation of the method 

The accuracy of the Riemann-sum approximation approach in (14) is validated by computing 
the steady-state solution for the velocity field. This is obtained by taking ^ = 0 in (5) and (6) which then 
reduces to the following ordinary differential equations. The implicit finite difference method has also been 
used to validate the Riemann-sum approximation approach, the advantage of this numerical procedure over 
others is that comparison can be made at both steady and transient state solution of the transport equations. 


^ + ±^_ M 2 l/ = -0 

CLR 2 R CLR 


(19) 


d 2 8 1 d8 _ q 

dR 2 RdR 


( 20 ) 


These are solved under the boundary conditions (8) to obtain the expressions for the steady-state 
velocity field, steady-state temperature field, steady-state skin frictions as well as the mass flow rate of 
the fluid. The solutions are respectively; 


I/O?) = C s / 0 (Mfi) + C 6 tf 0 (MR)+- 


e s (R) = i 


t 1s = ^| =M(C 5 / 1 (M)-C 6 /f 1 (M)) 

dR 'R=l 


( 21 ) 

( 22 ) 

(23) 


t As = - ^| R=A = M(C 6 /f 1 (AM) - C 5 I^AM)) (24) 

Qs = 2nf*RU(R,t)dR = 2n Q (A/, (AM) - I^M)) - J (A/s^AM) - K^M)) + ^ (A 2 - 1)] (25) 

The constants C 5 and C 6 in (25) are stated by: 


8 = 
El 
E2 
E3 


Ci 

C 2 

C 3 

C 4 

C 5 

C 6 


Vm 2 + s 

= s[l 0 (AVsPr)K 0 (VsPr) - I 0 (VsPr)K 0 (AVsPr)] 
= s[s(Pr - 1) - M 2 ][I 0 (8)K 0 (A8) - I 0 (X8)K 0 (6)] 
= M 2 [I 0 (AM)K 0 (M) - I 0 (M)K q (AM)} 

_ [fr 0 (VsPr)-ff 0 (AVsPr)] 


_ [/ 0 (AVsPr)-/ 0 (VsPr)] 
El 

[Jr 0 (Ag)-jf 0 cg)] 

E2 

U 0 a5)-/o(5)] 

E2 

{k 0 (Am)-k 0 (m)] 

E3 

[/ 0 (M)-t 0 (AM)] 


The numerical values of the velocity obtained using the Riemann-sum approximation approach, 
implicit finite difference method and those obtained from the exact solution of the steady-state choosing 
value of M — 2, A — 2. Is presented in Table 1. The comparison between the results, shows that at large time 
(steady-state) there is an excellent agreement between the Riemann-sum approximation approach and the 
implicit finite difference. 
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Table 1. Numerical values of the velocity obtained using Riemann-sum approximation approach, implicit 


finite difference and exact solution for different values of R 


t 

R 


Velocity 


Riemann-sum 

approximation 

Implicit finite 
Difference 

Exact 

solution 

0.2 

1.2 

0.0478 

0.0478 

0.0624 


1.4 

0.0652 

0.0651 

0.0870 


1.6 

0.0633 

0.0632 

0.0837 


1.8 

0.0434 

0.0434 

0.0553 

0.4 

1.2 

0.0608 

0.0608 

0.0624 


1.4 

0.0846 

0.0846 

0.0870 


1.6 

0.0814 

0.0814 

0.0837 


1.8 

0.0540 

0.0540 

0.0553 

0.6 

1.2 

0.0623 

0.0623 

0.0624 


1.4 

0.0868 

0.0868 

0.0870 


1.6 

0.0835 

0.0835 

0.0837 


1.8 

0.0552 

0.0552 

0.0553 

Steady state 

1.2 

0.0624 

0.0624 

0.0624 


1.4 

0.0870 

0.0870 

0.0870 


1.6 

0.0837 

0.0837 

0.0837 


1.8 

0.0553 

0.0553 

0.0553 


3. RESULTS AND DISCUSSION 

In order to have a clear insight of the physical problem under consideration, a numerical 
computation is performed using the mathematical laboratory software (MATLAB) to compute and generate 
graphs for the velocity field, temperature field, skin-frictions and mass flux for different values of 
the governing parameters, so as to comment on their relative contribution to the flow formation. In this work, 
two different cases of fluid are been examined these include air with Pr — 0.71 and water with Pr — 7.0. 
The effect of variation of the governing parameters M, Pr and X on the flow formations are presented in 
Figures 2-9. Unless otherwise stated, the value t = 0.2, M — 2 and X — 2 are selected arbitrarily to study 
the effect of various parameters on the flow behavior. The influence of Prandtl number and time on 
the temperature profiles is shown in Figure 2. It is revealed in Figure 2 that fluid temperature increases as Pr 
and t increases. It is concluded from Figure 2 that the reduction in fluid temperature is directly proportional 
to the decrease in thermal diffusivity. Figure 3 shows that an increase in the values of the Hartmann number 
M causes retardation to the fluid flow indicating the fact that the imposition of magnetic field slow down the flow. 


'W 


Pr = 0.71 
Pr = 7.0 


, t = 0.4, 0.6, 0.8, 1.0 


rm* * •jrrrsT 

k /•*•••. • •.....»•••••* # //I 




1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 

R 


0.08 

0.07 

0.06 

0.05 

— 0.04 

7D 0.03 
> 

0.02 

0.01 

0 

- 0.01 


--.. 

- 

.... 

—•— Pr = 0.71 

Pr = 7.0 - 


... 

• * ,* . • -. • • _ 

..— 

\ 


X*v\s 

- # \ 

% 

- - \ ' 

f M = 1.0, 1.5, 2.0, 2.5, 3.0 



1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

R 


Figure 2. Temperature distribution for different Figure 3. Velocity distribution for different values 

values of t (M = 2.0) of M (t=0.2) 


This remark is consistent with the physical fact that the Lorentz force that appears due to 
the interaction of the magnetic field and the fluid velocity resists the corresponding fluid flow, resulting in 
the velocity to decrease gradually. Figure 3 also indicates how the velocity field is affected corresponding to 
an increase in the values of Pr. We recall that an increase in Pr in Figure 2 signifies a fall in thermal 
diffusivity for the model under consideration. It is learnt from Figure 3 that when the thermal diffusivity of 
the fluid is reduced, the flow gets decelerated largely which may be attributed to the fact that a low thermal 
diffusivity leads to a corresponding decrease in the kinetic energy of the molecules of the fluid, which in turn 
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affects the fluid velocity adversely. Figure 3 also revealed that the fluid velocity increases with increase in 
time till it attains steady state. It is worthy to note that Pr has no effect on the velocity and temperature 
profiles at steady state. 

Figures 4 and 6 present variation of skin-friction at outer surface of the inner cylinder (R — 1) 
and the inner surface of the outer cylinder (R — A) respectively for different values of the Hartmann 
number (M). It is obvious from these Figures 4 and 6 that skin-friction decreases with increase in Hartmann 
number (M) on both surfaces for both cases of Pr. Fluid with Pr = 0.71 is observed to induces higher 
friction at both surfaces of the cylinders in comparison with fluid with Pr — 7.0 which is physically true, 
since higher velocity results in friction at the walls. In addition, it is clear that skin friction at both cylinders 
for air attains steady state faster than water. This suggests that if one considers to reduce the friction at the 
surfaces, fluids with higher Prandtl numbers like water (Pr — 7) should be considered. Variation of skin- 
friction profiles at the outer surface of the inner cylinder and the inner surface of the outer cylinder for 
different values of radii ratio (A) are shown respectively in Figures 5 and 7. It is obvious from both Figures 5 
and 7 that skin-friction on both cylinder increases with (A). A keen scrutiny of the Figures 5 and 7 reveal that 
the skin-friction is independent of time for Pr — 0.71 except as t tends to zero. 



Figure 4. Variation of skin friction (r^ for different Figure 5. Variation of skin friction(r-,) for different 
values of M values of A(M — 2) 



Figure 6. Variation of skin friction( t^) for different Figure 7. Variation of skin friction(r^) for different 

values ofM(A = 2) values of A(M = 2) 


Figure 8 reveals the mass flow rate plotted against Hartmann number (M) for different values of 
time (t). It is observed that the mass flow rate decreases as Hartmann number (M) increases but increases 
with time (t). Figure 9 illustrations variation of mass flow rate for different values of radii ratio (A). It is 
evident that mass flow rate increases with increase in radii ratio (A) for both cases of (Pr). It is worthy to 
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note that mass flow rate is constant for different values of radii ratio (2) in the case Pr = 0.71, but increases 
with time (t) for Pr — 7.0. Figures 8 and 9 lead us to conclude that the parameters Pr and M have significant 
contributions in regulating the amount of total discharge of fluid through the annulus and they may be 
suitably chosen to control the mass flux. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

t 



Figure 8. Variation of mass flow rate (Q) for 
different values of M(2 = 2) 


Figure 9. Variation of mass flow rate ( Q ) for 
different values of 2(M = 2) 


4. CONCLUSION 

A semi-analytical study is conducted to examine the role of magnetic field on an incompressible 
and an electrically conducting fluid filled within two coaxial cylinders. The Laplace transform technique 
and Riemann-sum approximation method have been used to obtain the solution of the governing equations. 
The influence of Hartmann number (M), Prandtl (Pr), radii ratio (2), and time (t) on the velocity field, 
temperature distribution, skin frictions and mass flow rate have been extensively discussed. The main 
findings in the present research are: 

a. It is found that an increase in Hartmann number (M) has a retarding effect on the velocity field, mass 
flow rate and skin-friction at both surfaces. 

b. It is worthy to conclude that an increase in the radii ratio (2) increases the skin-friction at both surfaces. 

c. Air is established to have higher fluid velocity, mass flow rate and skin-friction on both surfaces in 
comparison with water. 

d. Generally, Air (Pr — 0.71) attains steady state temperature faster, due to its higher thermal diffusivity 
(See Figure 3) in comparison with water (Pr = 7.0). 

e. Skin-friction at both surfaces and Mass flow rate are seen to be independent of time except at small value 
of time (t). 
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